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^ We show that the anomalously large top quark forward-backward asymmetry observed by 

^^ CDF and D0 can naturally be accommodated in models with flavor-violating couplings of a 

CN new massive spin-2 state to quarks. Regardless of its origin, the lowest-order couplings of a 

(— I spin-2 boson to fermions are analogous to the coupling of the graviton to energy/momentum, 

Ph 
I leading to strong sensitivity of the effects associated with its virtual exchange to the energy 

^ scales at hand. Precisely due to this fact, the observed dependence of the asymmetry on the tt 

invariant mass fits nicely into the proposed framework. In particular, we find a vast parameter 

^ space which can lead to the central value for the observed forward-backward asymmetry in the 
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OC high mass bin. 
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I. INTRODUCTION 

The Standard Model (SM) is a successful theoretical framework for describing ele- 
mentary particle interactions when confronted with experimental data. However, recent 
observations [H [2] by the CDF and D0 collaborations point to an anomalously large 
forward-backward asymmetry in ti production (A^^), significantly exceeding the SM 
prediction (see Ref. [3] for a review). 

Perhaps the most intriguing result is CDF's report [1] of a rise in A^p^ with the invariant 
mass of the ti pair (Mtf), 

AL = A%{Mti < 450 GeV) = (-11.6 ± 15.3)% , 
<gh = A%{Mu > 450 GeV) = (47.5 ± 11.4)% . 

In particular, A^-^ is almost 3cr away from SM prediction (all the relevant measurements, 



as well as the corresponding SM predictions are collected in Table III). 

This discrepancy invites a new physics (NP) explanation, and many models have been 
proposed to address the anomalously large A^p^. These models generally involve intro- 
ducing new scalar [5HT3] or vector [T3H2D] particles contributing to the ti production 
cross section in the s- and/or t-channel. While most of these models can easily raise the 
theoretical prediction for A^p^ to within la of the CDF measurement, it has proven to be 
extremely hard to address the central value of A^-^ while being consistent with existing 
experimental constraints. 

In this work we propose to explore another class of models involving new tensor (spin- 
2) particles around the weak scale with flavor-violating couplings to quarks. A simple 
effective field-theoretic (EFT) analysis reveals that the most general lowest-order cou- 
plings of a spin-2 state with quark bilinears are rather similar to the general-relativistic 
couplings of the graviton to energy /momentum. This leads to a strong energy dependence 
of the effects of virtual exchange of such states, which nicely agrees with the CDF ob- 
servations. In particular, we show that this framework can accommodate all of the CDF 
measurements over a wide range of parameter space while being consistent with existing 
experimental bounds. 

We treat the massive spin-2 particle as a low-energy signature of some unspecified UV 
physics. Among other possibilities, a low-energy effective theory of this type could arise 



from a theory of modified gravity [2T1426] or could describe a spin-2 resonance of a strongly 
interacting sector not far above the weak scale [271 - I29] . 

The rest of the paper is organized as follows: in Sec. [IT] we describe the relevant cou- 
plings of the spin-2 state to quarks. The existing experimental constraints are examined 



in Sec. Ill We study parameter space for tt phenomenology and discuss our results in 



Sec. IV Finally, we conclude in Sec. ^ 



II. LOW-ENERGY EFFECTIVE THEORY 

Studies of low-energy effective field theories for massive spin-2 particles are motivated 
in various contexts. On one hand, such theories can be useful for describing spin-2 QCD 
resonances, such as glueballs, which become long-lived in the limit of large number of 
colors [29] . On the other hand, massive spin-2 states frequently occur in the context of 
modified gravity. An incomplete list of examples from the latter category includes models 
with KK tower of gravitons, such as theories with large [211 122] or warped [231 El] extra 
dimensions, as well as the recently discovered class of purely four-dimensional, ghost-free 
models of massive gravity [25] [26]. Ratios of branching ratios to photons and to jets may 
be used to distinguish between the various possible underlying UV theories [30] . 

Regardless of the details of the UV theory, any consistent action for a complex, sym- 
metric spin-2 field /i^^ with mass M should reduce to the Fierz-Pauli [31] form at the 
linearized level 

^FP = -hi (D + M') h^'' + \ h^J (D + M') K - h\,d^d^h^^ + hl^d^d'^h^^ + h.c. (1) 

Furthermore, if the field is of gravitational origin, its couplings to matter are usually 
constrained by the Equivalence Principle to be universal. However, RS-type models with 
the SM fields localized differently along the bulk are an important exception to this rule. 
In the present work, we will be mostly interested in the implications of a massive, 
complex, spin-2 boson for the top quark forward-backward asymmetry. Resorting to an 
EFT approach, we will not make any assumptions about the precise origin of the spin-2 
field. Among other possibilities, h^^ could describe a bound state of some strongly coupled 
sector not far above the electroweak scale, or a non-universally interacting gravitational 
KK mode in some RS-like theory with complicated localization of matter. We will not 



attempt to construct any explicit model along these lines. 

At low energies, the sector of the theory describing interactions of h^^ with the quarks 
consists of operators of various dimension, suppressed by powers of some high scale, 
denoted by / below. At the zero derivative order, there is a single coupling of a spin-2 
field with a general quark bilinear, 

£4 3 Aij/i^^gigj + h.c, (2) 

where A is a coupling constant, {^, j} refer to quark flavor and the possible chirality of 
qi has been suppressed for simplicity; the quark fields correspond to the mass eigenstates 
after electroweak symmetry breaking. This interaction is similar in form to an ordinary 
Yukawa interaction of a SM singlet scalar; the only difference is in the spin-2 nature of 
the correlator 

For the purposes of studying the ti forward-backward asymmetry however, this operator 
can be expected to lead to effects similar to those of a color-singlet scalar exchange. 
The consequences of the latter for the ti forward-backward asymmetry have been studied 
extensively, see for example [5], with the conclusion that it can not generate a large 
enough asymmetry. Moreover, such inter-generation couplings involving the top quark 
are constrained not to be too large from the same sign top production cross section at 



the LHC (see Section III); we will thus ignore these operators below^. 

At the one derivative order, the most general couplings of /i^i, to a fermion bilinear are 
given by the following expressions, 

/:5 3-^V('5^' + ^'''r) + h.c., (4) 

Sfiu = jO^ijQLi {ifidy + 7^.(9^) qij + -&j^- {df^qLilv + d^qLil^i) qij + (L ^ R) , (5) 

T = XfjqLi^qLj + {L^R), (6) 

with arbitrary complex coefficients {a^'^, 6^'^, A^'"^}. For external (on-shell) fermions, the 
interactions in ^ can effectively be reduced via the Dirac equation to the non-derivative 



^ If hfj^i, is of gravitational/extra dimensional origin, / represents the quantum gravity scale. Then the 
coupling of h^i, to energy-momentum tensor leads to a natural suppression of the Yukawa coupling 
constant. In particular, A ^ rn/ f, where m is of the order of the mass of a heavier fermion present in 
the interaction. 



Yukawa couplings given in ([2]). Again, since we are interested in the leading order spin-2 
effects on fermion scattering, we will ignore the operators in ^ for the remainder of the 
paper. 

A priori, there are no constraints on the couplings a^'^ and b^'^. However, for a theory 
in which the spin-2 field is gravitational in nature, S^, should be related to the energy- 
momentum tensor. The fact that /i^jy interactions are non-universal, as well as couple 
different generations with each other does not rule out the possible gravitational interpre- 
tation of the theory - this can be accommodated e.g. in the framework of RS models with 
complicated matter localization along the bulk. In such cases the KK gravitons couple to 
quark flavor eigenstates in a diagonal, albeit non-universal way. This, upon rotation to 
the mass basis, results in the following constraints on the couplings 

^ij = -bij = Qij . (7) 

Although we remain completely agnostic about the origin of the spin-2 state, we will take 
these relations to hold in the analysis to follow; as we show below, restricting the parameter 
space in this way is already enough for generating the needed amount of asymmetry 
without running into conflict with other experimental bounds. Further constraints on the 
couplings come from the requirement of the invariance of the theory under the SM gauge 
group. We will not dwell on making the symmetry manifest, but will keep in mind that 
it implies some additional relations between the coupling constants. 

The key observation to make at this point is that even if the spin-2 state is not asso- 
ciated with any gravitational dynamics, its couplings to fermions are quite similar to the 
coupling of the graviton to the energy/momentum. This leads to a large sensitivity of the 
effects associated with the virtual exchange of these states to the energy scales at hand. 
In particular, we will find that this fact fits with the observed pattern of an increase in 
A^PQ with ti invariant mass. 

The range of validity of the low-energy effective theory is an important issue, since it 
determines the maximum energy to which the analysis performed below can be extrapo- 
lated. Given the complete ignorance of the UV completion of low-energy (massive) spin-2 
theories, the best we can do is to make an educated guess of the relevant energy scale. 
For massless general relativity (GR) in 4D, the analysis of graviton loop corrections yields 



the following expansion parameter |32l 133] , 

where E represents the typical energy scale of a process under consideration, This is in 
complete analogy to what one finds for low-energy nonlinear sigma models, once Mpi is 
replaced by the pion decay constant, f-„. 

The massive spin-2 representation of the Poincare group propagates three additional 
(one helicity-0 and two helicity-1) degrees of freedom on top of the two helicity-2 modes 
of the massless theory. In a general nonlinear completion of Fierz-Pauli action, the strong 
coupling of the helicity-0 mode is usually responsible for the cutoff of the effective theory 
- in complete analogy to massive non-abelian theories, where the strong coupling of the 
longitudinal PF-bosons leads to the necessity of UV completion at a scale ~ ATcmw/g (here 
rriyy and g refer to the W^-mass and the SU{N) coupling constant). Due to the higher 
spin structure, non-linear sigma models involving the longitudinal modes of massive spin- 
2 theories are usually different in nature than those for spin-1^. This leads to a UV cutoff 
which is in general sensitive to the nonlinear completion. For example, in theories of 
massive (four-dimensional) GR with a graviton potential, the cutoff usually comes out to 
be a certain geometrical mean of the scales M and /. A specific class of potentials which 
avoids the propagation of ghosts in the theory [2S| leads to a sigma model with higher 
UV cutoff, compared to theories with a more general potential [M]; some other possible 
UV/nonlinear completion (e.g. a completion beyong the potential, or one which relaxes 
the requirement of reproducing four- dimensional GR in the M — )■ limit) can therefore be 
expected to yield yet more different sigma models. Another possibility is that new physics 
regulating the low-energy theory in the UV kicks in somewhat below the strong coupling 
scale^, however it does not distort the spin-2 exchange effects up to higher energies. 



^ In particular, a longitudinal scalar mode of a massive spin-2 boson is more strongly coupled than that 
of a massive spin-1 particle. Mathematically this can be traced back to the piece in the massive spin-2 
propagator which grows fastest with momentum, 

f^fj. f^v f^a f^lS 



{hfiuhaji) D 



M'^{k^-M^)' 



^ In extra dimensional examples, one might imagine the higher KK modes softening loop effects, thus 
providing a completion at intermediate energies - up to the fundamental quantum gravity scale. 
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As already emphasized, below we will not be concerned with the nature of the un- 
derlying UV/nonlinear theory and will expect the low-energy description to be valid up 
to scales somewhat above the scale / = min{f /gi, M} where Qi collectively denotes all 
cupling constants in (|2| - (|6|. 

III. EXISTING EXPERIMENTAL CONSTRAINTS 

We now turn our attention to the analysis of the existing experimental constraints on 
the massive spin-2 model considered in this work. Note that the main motivation of the 
present work is an illustration of strong energy-dependence of the effects of virtual spin- 
2 exchange and its possible relevance to observations that exhibit these effects. In this 
section we are concerned with preliminary estimates of the bounds on the spin-2 parameter 
space due to various experimental constraints - just to show that the mechanism can be 
viable, or even robust. Of course, many additional studies need to be performed (e.g. a 
closer inspection of LHC bounds or studies of spin correlations of top quark daughters) 
for a complete phenomenological analysis. Wherever loop contributions are involved, we 
make the most conservative assumptions for their magnitude just to show that even with 
these overly restrictive assumptions, there still is a vast parameter space for addressing the 
ti forward-backward asymmetry. We leave a more detailed analysis of the experimental 
bounds for a future study. 

Below we determine bounds from LEP, electroweak precision data, same-sign top-quark 
pair production, Bd — Bd mixing and dijet production at the Tevatron. As can be seen 
below, the effects of these constraints on the model are mild. We note however, that a 
recent analysis |35] of the kinematic and dynamical aspects of the relationship between 
the asymmetries A^p^ and A^p^ measured by D0 favors new physics (NP) models that 
produce more right-handed than left-handed top quarks. Even if the left-handed sector 
is taken to be suppressed, the right-handed sector can still accomodate the asymmetry, 
as we will show below. Such a suppression is further motivated by other experimental 
constraints, such as Bd — Bd mixing. It is also worth noting that the ti production cross 



section, Eqs. (13)-(16), of the model at hand is symmetric under the exchange of the left- 



and right-handed couplings. 




FIG. 1: One possible diagram for 4-jet production at LEP. The dashed hne represents 

the spin-2 particle exchange. 

A. LEP Constraints 

The LEP constraints depend on how the massive spin-2 state couples to electrons 
and final state quarks. Since we are only interested in generating a large A^p^ at the 
moment, we can take gut ~ C^(l) while allowing freedom for the couplings to leptons. In 
this scenario therefore, we do not anticipate any bounds from direct production at LEP. 
However, due to SU(2)w symmetry, g^^ ^i = (^f ^i where {m, t] stands for any combination 
of u and t, and a light spin-2 particle could lead to anomalous 4-jet events as shown in 
Fig. ^ For examining this bound, we can implement the results in the literature (see e.g. 
Sec. IV. A. of Ref. [17] and references therein) in our model. The amplitude for a 4-jet 
final state in the present model is suppressed by extra factors of E'^ j p compared to the 
case of a new scalar or vector field. Here E denotes a relevant energy scale in the process. 
The final operating energy of LEP II is 209 GeV and for the parameter space considered 
below / is the highest scale in the processes at hand. Even if we conservatively take these 
suppression factors to be 1 therefore, the bound on the mass of the spin-2 particle from 
LEP is quite mild, M > 100 GeV. 

B. Electroweak Precision Tests 

Electroweak precision data (EWPD) can provide strong constraints on models that 
attempt to explain the ti forward-backward asymmetry [36j . Corrections to EW precision 
observables due to the intermediate spin-2 state do not occur at tree level assuming it does 
not directly couple to the EW gauge bosons. At the one-loop level there is a contribution 
to the dimension-4 operator, CzuuZ^u'y'^u, arising from the diagram in Fig. [2j As shown 
in Ref. |37], the most stringent constraint on Czuu comes from atomic parity violation 




FIG. 2: The loop contributing to the EW precision observables such as Tz and Qw- 

experiments. The experimental and SM values for Qw in cesium atoms in the 2010 
PDG |38] can be turned into a bound on the NP contribution to the coefficient of the 
dimension-4 operator, IC^J^I < 1.3 ■ 10^^. An estimate of the spin-2 contribution to 
Qw{Cs), ignoring left-handed couplings, yields 

C. Single- Top Production 

Single-top, spin-2 production via the reaction ug ^ t h^^ is PDF enhanced at the LHC 
and PDF suppressed at the Tevatron relative to spin-2 mediated tt production, which has 
a qq initial state. The phenomenology of spin-2 production can be classified into two 
categories. In the ffist case h^^, is stable on collider time scales, which is predicted in 
large extra dimensions scenarios. The decay signature of this reaction - one b-tagged jet, 
one lepton, and high missing transverse energy (MET) - is not an event that is currently 
selected in single-top searches at the LHC. Single-top searches thus far always contain at 
least 2 jets or 2 leptons. The other scenario is that the spin-2 particle decays immediately 
upon production, which is the case in warped extra dimensions scenarios. In this case, 
bounds from single-top production can be avoided by making the branching ratio of the 



spin-2 particle into a uu pair small. As we will shown in Sec. IIIF the coupling guu, which 



controls the size of Br^h^u — > uu), is constrained to be small from dijet bounds. 



D. Same-Sign Top-Quark Pair Production 



The t-channel models aiming to explain the A^p^ asymmetry can be strongly con- 
strained by limits on the same sign top quark pair production at the LHC The ATLAS 
collaboration places limit on the production cross-section at au < 1-7 pb [53]. In our 
model, the production arises from processes shown in Fig. [3J These diagrams contribute 
to the coefficient of the effective 4-quark operator responsible for tt pair production, 
C{tR'~i^UR){tR'-)fj,UR). The CMS collaboration reports the bound on this coefficient of 
C < 2.7 TeV"' 



The tree-level diagram. Fig. 3a, constrains the combination of couplings dutidui)* / P ■ 
We have chosen to taken g!^^ = 1 in ([s]) so that the tree-level cross section yields a bound 
on gj^, which is given in Table ^ 

Since the spin-2 propagator contains pieces such as MHk^-A.p) ^ loops containing spin-2 
particles are highly divergent. The most conservative estimate (i.e. neglecting the possi- 



bility of some derivatives acting on the external fermions to reduce the energy-dependence 
of the finite part of the loop integral) of the one-loop contribution to same-sign top quark 
production yields 



C 



\9^ 



ut 



9l 



s s 



167r2 /4M8' ^^°^ 

where i = u, t and s is the partonic center of mass energy. Since the PDFs drop sig- 
nificantly for the momentum fraction greater than 0.3, we can estimate s = 10% of the 
LHC running energy, s = 700 GeV. For / = 350 GeV, M = 500 GeV and g^^ = 1, we get 
C ~ 3 ((7jf ) TeV~^. This leads to the bound gj^ < 0.9. Note that the bounds get weaker 
for a larger value of / and/or M. 

The bound from the non-derivative interactions in Eq. ^ can be easily obtained in a 
similar fashion. But since this coupling doesn't play much role in our analysis of A^p^, we 



M [GeV] 


100 


200 


300 


400 


500 


600 


700 


„R 

9tu 


0.05 


0.04 


0.03 


0.03 


0.03 


0.03 


0.02 



TABLE I: 



CL upper limit on l^f^l from a search for same-sign top-quark pair 
production by the ATLAS collaboration. I^f^^l is fixed to be 1, and / is set by requiring 



that ^high 



0.475. 
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(a) Tree-level diagram. (b) One of the possible 1-loop diagrams. 

FIG. 3: Feynman diagrams giving rise to the same sign top-quark pair production. 

will not dwell on the value of this bound. Put differently, we can take the coupling Xij to 
be negligible while still producing a large A^p^. 

Single-top, spin-2 production where the spin-2 particle immediately decays into tu 
or a ti pair can also contribute to same-sign top production. We aproximate (j{ug — )■ 
ttj) ^ cr{ug — > th^y) X [Br{h^u -^ tt) -\- Brih^y — )■ tu)) as we are only interested in 
producing a quick estimate of the cross section. The ATLAS collaboration [39j imposes 
a cut, \ri\ < 2.5, when selecting lepton and jet candidates assoicated with same-sign top 
production. We imposed this cut when calculating the spin-2 contribution to this cross 
section. We neglect Br{h^u -^ tu) since this branching ratio is proportional to Ifi'tuP, 
which is constrained to be small from same-sign top production in other channels. For 
M = 400 GeV and / = 1 TeV, the ATLAS bound on same-sign top production, au < 1-7 
pb, yields the constraint, \gut\'^Br{h^i, -^ tt) < 0.97. 

E. Bd - Bd Mixing 

Non-zero off-diagonal couplings can lead to flavor changing neutral currents (FCNCs). 
Here we focus on Bd — Bd mixing as the most restrictive bounds on g^^ are expected to 
come from this process. The spin-2 contributions to Bd mixing can be described by the 
following four-quark operator 



Qi= [dLl>.bL)[dLrbL). (11) 

In general, other operators contribute to Bd mixing as well. The coefficient of the operator 



in Eq. (11 ) is constrained to be smaller than (9(10 ^^) GeV [H]. These bounds constrain 



the couplings relevant for A^p^ in our model, in particular there is a constraint on g^^. 
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FIG. 4: A tree-level contribution to i?^ — B^ mixing. 



This is due to the fact that by SU{2)l symmetry, g^^^ = g^^. The contribution to this 
operator arises at tree-level from Fig. |4j A quick dimensional analysis estimate reveals 
that the contribution from such a diagram is 



9ut \9tu) 



mt 



^gutigtuY^O'^'GeV- 



P M2 

where we estimate rrib/M ~ 10^^ and take / = 500 GeV. Thus we see that for gutigfu)* ~ 
(9(0.01), the constraints from Bd — Bd mixing can easily be satisfied. As in the case of 
same sign top-quark pair-production, the constraint gets weaker for larger values of /. 
Similar constraints hold for the right-handed couplings. However, there is no symmetry 
relating g^^ to g^^. Thus there is more freedom available in the right-handed sector to 
address the top forward-backward asymmetry. 



F. Tevatron Dijet and Top Width Constraints 

As emphasized above, the derivative couplings of the spin-2 field to light quark pairs 
grow with energy, which can lead to strong constraints on the couplings at large invariant 
mass, Mjj. Following Ref. [17], we obtain the Tevatron dijet bounds from the CDF 95% 
CL upper limits on the product of an RS graviton (G*) production cross section x its 
branching ratio to dijets {B) x acceptance {A)\ see Table 1 of Ref. |12]. The results with 



M [GeV] 


300 


400 


500 


700 


900 


1100 


1300 


R 

iJuu 


0.08 


0.15 


0.16 


0.18 


0.17 


0.12 


0.09 



TABLE 11: 95% CL upper limit on l^f^^l from dijet constraints at the Tevatron. \g^f-\ is 
fixed to be 1, and / is set by requiring that v4"jgj^ = 0.475. 
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B-A = l are collected in Table [h} The CDF analysis uses the RS parameter k/Mpi = 0.1, 
which translates into M ^ 0.383/ in the notation of this paper (with the convention that 
the largest dimensionless coupling is set to one). The dijet cross section in the present 
model is therefore related to the RS graviton cross section by a^,^ = C^ac*, where 
C = \guu\ M/0.383f. To obtain a meaningful bound, / is taken to be of the value required 
to produce an asymmetry in the high mass bin of 47.5% with \g^f\ = 1. crj^^ is estimated 



to be the s-channel NP cross section, cr^^, see Eq. (15), includes additional terms due 
to a finite top-quark mass; Fi{x, 0) should be used be for light quark production cross 
sections as opposed to Fi{x, y). This is an overestimate of the actual dijet cross-section, 
which is produces a conservative estimate of the bound. 

g^^ does not play a role in generating a large asymmetry. However, g^^ along with 



TUU 



gj^ are important for increasing the total cross section from the SM value up to what is 



measured at the Tevatron (see Table III for the relevant measurements and the corre- 
sponding SM predictions). A smaller value of g^^ requires a larger value of g^ to produce 
the same cross section. A larger gf^ and a smaller g^^ could be expected in an RS model; 
localization of the top quark close to the IR brane leads to a large coupling, while the 
light quark couplings are relatively suppressed as they are localized in the bulk closer to 
the UV brane. 

The dijet invariant mass and angular distributions measured the LHC may very well be 
more constraining than the Tevatron measurements. However, we do not consider LHC 
dijet data because these measurements only constrain the coupling g^u in this model. 
As previously noted, guu plays no role in generating an asymmetry so constraining this 
coupling to be smaller does not affect the goal of this work. LHC dijet measurements will 
be important for constraining a more complete model, and we leave this work for a future 
publication. 



When M < rrit, the non-standard top quark decay F (t — > m h^,^) is allowed, 
r.NP_ \9ut + gtu\'ml f^ M'Y f M'\ 

The recent D0 measurement [13], F^ = (2.00]'^q;43) GeV, constrains how large \gut\ / f can 



be in the low mass region; see Figs. 8a and 11 for bounds on the spin-2 parameter space 
from Ff. 
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FIG. 5: Most important spin-2 contributions to tt production at the Tevatron. 
IV. TOP QUARK FORWARD-BACKWARD ASYMMETRY 
A. Calculation of the Cross Section 

In the present section we study the effects of an intermediate massive spin-2 state 
on the top quark forward-backward asymmetry at the Tevatron. Calculations of the 
virtual exchange of a massive spin-2 graviton in the context of theories with large extra 
dimensions were done by Giudice et al. in [33j. However, the results of Ref. [33] should 
be applied to collider phenomenology with care as they are only valid in the limit (in our 
notation) p ~ M^ ^ s ^ m\. Here we extend the results of |33| by assuming there is 
no hierarchy between the aforementioned scales. 

The weighted average (the average over initial spins and colors and the sum over final 
spins and colors) of the amplitude squared for the qq — )■ tt scattering from a spin-2 
t- channel exchange is 



(IMI^) ^' 



128/4 ('(t-M2)Vr2M2 
C| / m^ m^ m^ m^ 



.^' [^' + 5^^^ + S^^v 



(13) 



where s, i are the Mandelstam variables in the parton center-of-momentum frame, q = 
{u, c} is assumed to be massless, and F denotes the width of the spin-2 resonance. The 
functions Fi are polynomials in x = i/s and y = m'f/s, and are defined in Appendix [AJ 
The Cj's are combinations of couplings and are also given in Appendix [A} 

The interference of the t-channel spin-2 exchange with the SM leading order (LO) 
gluon exchange gives 



27ra,s2 (t - M2) / mf mf \ 

(2 Re (MMsm)) = . / ^ . ^ V ^2 r^ + -W2^'' + -Wi^^' 



(14) 
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while the weighted average of the amphtude squared for the qq — )■ ti scattering from a 
spin-2 s-channel exchange is 

(|A^/) = -, o r {C^ Fi2 + Ca Fi3 + Cs F14) , (15) 

^' '^ 128/4 ((s - M2)2 + r2M2) ^ ' ^' ^ ' 

where q = {u, d, s, c, b}. The exchange of any color-singlet particle in the s-channel 
can not interfere with color-octet gluon exchange in the SM. However, s-channel spin- 
2 exchange can interfere with the exchange of a spin-2 particle in the t-channel with 
q = {u, c}, 

(2 Re {MMD) = —^ -^ -y- -^ ^ (16) 

1152/4 ((s - M^f + r2M2) Ut - MY + T^mA 

X _^6 [Fi, + —^F,, + —F,,j + Cj ^Fi8 + — Fi9 + — F20 
Our results are consistent with what was found in Ref. 



B. Tevatron Measurements and SM Predictions 

We start out by reviewing the recent observations of the anomalously large top quark 
forward-backward asymmetry, A^pg. The experimental evidence for contributions to A^p^ 
from physics beyond the SM is as follows. The CDF collaboration measured pj the 
asymmetry to be (20.0 ± 7.0)%. A recent D0 analysis |i2j yielded the value A"^ = 
19.6]']g'5%, in good agreement with the CDF measurement. D0 also reports a forward- 
backward asymmetry based on the rapidity of the leptons from top quark decays of A^p^ = 
(15.2 ± 4.0)% compared with the small SM value (2.1 ± 0.1)% calculated using MCONLO. 
All uncertainties have been added in quadrature. In addition, the CDF collaboration 
reports f3J that the asymmetry rises with the invariant mass of the ti system, with A^-^ = 
A§s {Mti > 450 GeV) = (47.5 ± 11.4)%, and Ag^ = Af^ {Ma < 450 GeV) = -(11.6 ± 
15.3)%. 

Despite recent improvements in the SM calculations, the asymmetry in the high mass 
bin is still close to three standard deviations away from the SM value. The central value 



■^ For example, ^12(2:, 0) — G4{x) and F4{x, 0)/18 + ^15(2;, 0)/6 ~ Gii{x), where Gi{x) are given in the 
appendix of [33]. 
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of a next-to-leading order plus next-to-next-to-leading logarithm (NLO+NNLL) QCD 
calculation of A^p^ is 7.3to/i% |16]. Recently calculated electroweak Sudakov (EWS) 
corrections enhance the QCD asymmetry by a factor of 1.041 to 7.7% [H], while fixed 
order electroweak contributions add an additional 1.6% to the asymmetry [ITj. The 
overlap between the EWS corrections and the fixed order EW contributions is estimated 
in Ref. [S] to be ~ 0.5%. This yields a total SM prediction of Aff = 9.2+|^ %. Similarly, 
combining the QCD predictions of Ref. [M] (which use MSTW2008 PDFs |19]) and EW 
effects calculated in Refs. [HI HT], the total SM prediction for A|.^ in the low and high 
mass bins is 5.4to'Q % and 14.1^20 % respectively. The 2.8a deviation from the SM in the 
high mass bin may be taken as a signal of new physics (NP). 

The total cross section for ti production was recently measured [35] by CDF to be 
att = (8.5 ± 0.9) pb. This measurement is consistent with the value reported [50] by D0, 
att = 7.78toll pb. Cacciari et al. [51] calculated the total cross section at approximate 
NNLO QCD to be 6.722l|J;4i^ pb. The EWS correction factor for this observable is 
7^i = 0.98 [M]. For Mti < 450 GeV the NLO+NNLL SM prediction for a^*^ is 4.23 
pb, while for Mti > 450 GeV a^^^ = 2.40 pb [17]. TZt in these bins is 0.985 and 0.973 
respectively^. The measured and predicted values for ti observables at the Tevatron are 
summarized in Table I 



Observable 


Measurement 


SM prediction j44j 


Att 

■^FB 


(20.0 ±4.7)% [E! 


9.3l|^% 


Att 

^high 


(47.5 ±11.4)% m 


i4.ii;^% 


Att 


-(11.6 ±15.3)% [H 


5.4l°j% 


c^tt 


(8.5 ±0.9) pb [15] 


6.59lJ]i^ pb 



TABLE III: Measurements and predictions for observables in tt production at the 

Tevatron. 



^ We thank Mike Trott for the computation of the EWS correction factors for ati in the low and high 
mass bins. 
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C. New Physics Results 

Here we show that tree-level exchanges of a massive spin-2 particle can contribute 
significantly to v4"^, ameliorating the tension between measurements at the Tevatron 
and SM predictions. Following Ref. [TTj, we define a partonic level asymmetry, 

NP _ NP SM 

aNP+SM _ ^F ^B I aSM ^ /-lyx 

"^FB - (^5A/)^^ + ^NP + "^FB ^SM + ^JVP ' ^-^^^ 

which is to be compared against the binned partonic asymmetries reported in [3]. For 
later convenience, we define 

NP _ NP 
aNP _ ^F ^B /-I ON 

^FB - (^5M)^^ + ^NP i^^) 

as the NP contribution to the top-quark forward-backward asymmetry. We use state-of- 
the-art predictions for the SM quantities Af,^ and o"^^^, and LO predictions for the NP 
corrections. The partonic NP cross sections are convoluted into hadronic cross sections 
using NLO MSTW2008 parton distribution functions (PDFs) [49j. The factorization and 
renormalization scales are taken to he fi = rrit = 173.1 GeV, while Og is set by the MSTW 
fit value: as (Mz) = 0.12018. The width of the spin-2 state is taken to be a tenth of its 
mass, r = M/10. We take all of the couplings to be real, and ignore contributions from 
Fi4:{x,y) as it is numerically negligible. 

The total asymmetry as a function of / is shown in FigJ6 Here U^J = 1 and we have 
set all other couplings to zero. The thick, dashed line is the central value as measured 
at the Tevatron with la error bands. The dotted line is the SM prediction. From left to 
right (at the top of the plot), M decreases from 700 GeV to 100 GeV in steps of 100 GeV. 
For intermediate values of /, the destructive interference with the SM exceeds the pure 
NP contribution, decreasing the asymmetry. As expected, the NP decouples for large /. 

Fig. [7] shows the effects of NP on ati as a function of /. The cross section for t-channel 



7b 



NP is shown in Fig. 7a with u^^J = 1 and all the other couplings set to zero. While Fig. 
shows the effect of s-channel NP as a function of / with the only non-zero coupings being 
Iduul — Idttl — 1- I^ Fig. 7a, M again monotonically decreases from 700 GeV to 100 GeV 



from left to right in steps of 100 GeV. In Figure 7b however, the ordering is not as simple, 
M = {100, 200, 300, 400, 700, 500, 600} from left to right. 

A global fit of the spin-2 model to the CDF measurements ^high; ^low cindcrti was per- 
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800 

/ [GcV] 

FIG. 6: A^pQ vs. /. \gut\ = 1 cind all other couplings are zero. The thick, dashed line is 

the central value as measured at the Tevatron with la error bands. The dotted line is 

the SM prediction. From left to right, M decreases from 700 GeV to 100 GeV in steps of 

100 GeV. 



formed using the method of least squares assuming the measurements are uncorrelated^. 
The scale / was fixed to be 1 TeV and g^^, g^^, and (7^ were left as free parameters for 
a given M. The results of the fit are shown in Fig. |8} The 1 and 2a confidence regions 
of allowed parameter space are shown in green and yellow respectively. The black line 
corresponds to A^-^ = 47.5%. This is not necessarily the best-fit value. Experimentally 
disallowed parameter space due to the constraints from same-sign top production at the 
LHC, EWPD, and the width of the top quark are shown in blue, red, and brown respec- 



tively (see Section III for a detailed discussion of experimental constraints on the model). 
The spin-2 model is able to hit the central value of the forward-backward asymmetry in 
the high mass bin in a large region of the parameter space. 

Fig. [9a] shows the binned asymmetry predicted by the spin-2 model for M = 350 
GeV, \g^^\ If = 2.36/TeV, and all other couplings set to zero. The CDF measurements 
with error bars are also shown. The purple band represents the theoretical uncertainty 
from varying the factorization scale in the range n = {mi/2, 2mt}. This combination of 



^ Of course, the measurements actually are correlated, but this should not affect the conclusions we are 
able to draw from the fit in any qualitative way. 




800 inoo 

/[GeVl 



loo 



(a) i-channel NP 



FIG. 7: att vs. /. t-channel NP is shown in Fig. 



/ |(i<'Vl 

(b) s-channel NP 



7a 



with \g^ 



ut 



1 and all other 



couplings set to zero. Fig. [7b] shows s-channel NP with the only non-zero coupings being 
1. In Fig. 



Huu 



i/tt 



7a 



M decreases from 700 GeV to 100 GeV from left to right in 



steps of 100 GeV. In Figure 7b however, the ordering is not as simple, 
M = {100, 200, 300, 400, 700, 500, 600} from left to right. 



parameters hits the central value of A^^ in the high bin and is within la of the central 
value in the low bin. The sum of the SM LO prediction plus the contribution from the spin- 



2 model with the same parameters as those used in Fig. 9a for the binned differential cross 



section, dau/dMu^ is shown in Fig. 9b Again, the purple band represents the uncertainty 

are 



in the PDF factorization scale, and the CDF measurements, as reported in Ref. 
also shown. The high bin values do not agree with the CDF measurements. However, 
we have not taken into account any detector acceptance effects. The deconvolution to 
the parton level done by CDF assumes the SM. As shown in Ref. [1^, model-dependent 
acceptance effects can reduce the cross section by as much as a factor of ~ 1/2 in the 
high bins. Furthermore, the total cross section reported in Ref. [22] is (Tu = 6.9 ± 1.0 pb, 
which is lower than the most recent measurements from both the CDF [IS] and D0 [HD] 
collaborations. It is reasonable to assume that detector acceptance effects and the known 
increase in the measured value for au could account for the disagreement in the high mass 
bins for dau/dMu- 
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FIG. 8: Results of a global fit of the spin-2 model to Tevatron observables. 

^high ~ 47.5% is shown in black. The 1 and 2a confidence regions of allowed parameters 

are shown in green and yellow respectively. The blue, red, and brown regions are 

disfavored by constraints from same-sign top, EWPD, and the width of the top 

respectively. 
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FIG. 9: Prediction from the spin-2 model for A"^ and dau/dMu with M = 350 GeV. 

The purple band represents the theoretical uncertainty from varying the factorization 

scale in the range /i = {mt/2, 2mt}. This example hits the central value of A^p^ in the 

high bin and is within la of the central value in the low bin. Detector acceptance effects 

and the known increase in the measured value for ati could account for the disagreement 

in the high mass bins for dati/dMti- 

D. LHC Measurements, Predictions, and Results 



There is no forward-backward asymmetry at the LHC because of its symmetric initial 
state, pp, as opposed to the pp initial state at the Tevatron. However, the same underlying 
physics that leads to ApB results in a charge asymmetry at the LHC, which we define as 

a (Ay2 > 0) - a (A|/2 < o) 



A^r 



c 



a (A|/2 > 0) + a (Ay2 < o) ' 



(19) 



where Ay^ is the difference of the squares of rapidities of the top quark and anti-top 
quark, Ay^ = yf — y^. The CMS collaboration reports [53] the charge asymmetry of 
Al = (-1.3^4:2) %> which is consistent with the SM prediction A^ = (1.15 ± 0.06) 
The value of A^ reported [55] by the ATLAS collaboration, Af, = 



-2.4 ±2.8)%, is 



consistent with the measurement of CMS. 

ATLAS recently measured [56] the top quark production cross section with ^/s = 7 TeV 
to be att = 176!^j;4 pbj which is consistent with the CMS observation au = (154 ± 18) 
pb [5?J. A QCD prediction |46j at approximate NNLO (using IPIscet) yielded au = 
155^12 pb- The EWS correction factor [S], TZt = 0.98, is used to compute the full SM 
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20 30 40 

A~^(i\/„->450GeV)[%] 

FIG. 10: NP contributions to the inclusive charge asymmetry at the LHC and the 

forward-backward asymmetry at the Tevatron in the high mass bin. The pink band is 

the la error of the measurement ^wh' ^^"^ ^^^ dashed hne is the difference between the 

measured value of A^jgh and the corresponding SM prediction. The blue and red curves 

are predictions of the spin-2 model for M = 100, 200 GeV respectively. 



prediction, a^ = 152;;'^}2 pb- 

shows A^ at the LHC as a function of A'^p^ at the Tevatron for the spin-2 



Fig. 



10 



model at hand. The pink band is the la error of the measurement A^-^, with the central 
value given by the dashed line. The blue and red curves are predictions of the spin-2 
model for M = 100, 200 GeV respectively. These predictions are not within la of both 
measurements simultaneously. This is a generic feature of any model that attempts to 
explain Af,^ [581 ES] • 

If the only NP is a single spin-2 field, then the cutoff of the effective theory should 
be at least as large as the center-of-mass energy of the experiment. The most optimistic 
estimate of the cutoff is A ?^ 47r/. As argued in section |ll|, the cutoff is likely to be 
smaller than this estimate. Fixing / to be 1 TeV, Qut should be less than <, 1.8 if the 
effective theory is to be valid up to 7 TeV, as opposed to gut ~ 6.4 for a/s = 1.96 TeV 
at the Tevatron. Adding heavier fields to the effective theory could raise the cutoff, as 
well as lead to a qualitatively different relation between the forward-backward and charge 
asymmetries measured at the Tevatron and LHC. Preliminary analysis has shown that 
interference between the virtual exchange of a lighter spin-2 particle and a heavier spin-2 
state could ameliorate the tension between the measurements of A^ and A'^^^. However, 
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FIG. 11: Fit to y4**gh (CDF) and a„- (ATLAS). / = 1 TeV and g^^ is a free parameter for 

a given M. The 1 and 2a confidence regions of allowed parameter space are shown in 

green and yellow respectively. The black line corresponds to A^igh — 47.5%. 

Experimentally disallowed parameter space due to constraints the width of the top 

quark is shown in brown. 

a more comprehensive study of this effect is needed for quantitative predictions, and we 
leave this for a future study. 

A global fit of the spin-2 model to the CDF measurement ^wh ^'^d the ATLAS mea- 
surement Gti was performed using the method of least squares. The scale / was fixed 
to be 1 TeV and g^^ was left as a free parameter for a given M . Again, g^t should be 
less than <, 1.8 if effective theory is to be valid up to 7 TeV. The results of the fit are 



shown in Fig. 11 The 1 and 2a confidence regions of allowed parameter space are shown 



in green and yellow respectively. The black line corresponds to A 



high 



47.5%. This 



is not necessarily the best-fit value. Experimentally disallowed parameter space due to 
constraints the width of the top quark is shown in brown. As was the case with the fit 
to Tevatron observables, the spin-2 model is again able to hit the central value of the 
forward-backward asymmetry in the high mass bin in a large region of the parameter 
space. 

There is an additional contribution to the tt production cross section at the LHC from 
single-top, spin-2 production where the spin-2 particle immediately decays into wt. Again, 
we approximate the 2 — )■ 3 cross section as o{ug — )■ ttu) ~ o{ug — )■ th^^) x Br(h^y — >• 
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ut). We used the ATLAS collaboration's [56j cut, |?7| < 2.5 when selecting muon and 
jet candidates associated with ti production, when we calculated the spin-2 contribution 
to this cross section. For M = 200 GeV and / = 1 TeV, requiring ati to be within 
la of the measured value limits \gut\ to be less than 1.55 assuming Br{h^i, — )■ ui) = 1. 
Alternatively, one may allow \gut\ to reach its maximum allowed value in the effective 
theory of approximately 1.8 by requiring Br{h^y — )■ ut) ^ 0.9 The contribution to au from 
this channel falls with the mass of the spin-2 particle such that theoretical considerations 
quickly become the dominant constraint on \gut\- 

E. Comments on Differential Measurements 

At the time this paper was submitted for publication, new measurements of the charge 
asymmetry at the LHC at the differential level were reported by the ATLAS [SS] and 
CMS \EU\ collaborations. CMS also recently measured |ni] the normalized, differential 
ti production cross section. By normalizing the differential cross section to the total 
cross section, certain systematic uncertainties and all normalization uncertainties cancel 
out, leading to a particularly precise measurement. In addition, the CDF collaboration 
updated [62] its analysis of the forward-backward asymmetry at the Tevatron to include 
the full Run II dataset. It is observed that A"^ has an approximately linear dependence 
on both Mft and Ay. In this section, we discuss the spin-2 model's predictions for these 
differential measurements. 



Fig. [12] shows predictions of the spin-2 model for the forward-backward asymmetry at 
the Tevatron, the normalized differential top quark production cross section at the LHC, 
and the charge-asymmetry at the LHC. The blue, red, and green lines correspond to a 
spin-2 mass of {100, 200, 300} GeV and a coupling Qut/f = {0.85, 1.38, 1.84} TeV~\ The 
dashed lines are the SM values for these observables. These calculations were made using 
FeynRules [63] interfaced with MadGraph 5 [6^ . 

As was the case for the inclusive charge asymmetry, these predictions are not simul- 
taneously within la of both the differential charge asymmetry and the forward-backward 
asymmetry measurements. To the best of our knowledge, this is a generic feature of any 
model for A"^ proposed before this paper was submitted for publication, see for example 
Refs. [HHl |59]. After this paper was submitted for publication, Drobnak et al. discov- 
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FIG. 12: Predictions for forward-backward asymmetry at the Tevatron (upper left), 

normalized, differential cross section at the LHC (upper right), and charge asymmetry 

at the LHC (bottom). The blue, red, and green lines correspond to a spin-2 mass and 

coupling, gut/f, of (100 GeV, 0.85 TeV-^), (200 GeV, 1.38 TeV"^), and (300 GeV, 1.84 

TeV~^) respectively. 



ered [HS] a class of models that can accommodate both measurements simultaneously. 
Based on the results of [65] , a spin-2 field that is charged under certain representations of 
the SM gauge group may produce a large A^^ and a negligible A". Verifying this claim 
is beyond the scope of this work, but it will be investigated in a future project. 

The shape of the normalized, differential cross section predicted by the spin-2 model 
does not agree with the CMS measurement. However, as was the case for the differential 
cross section at the Tevatron, we have not taken into account any detector acceptance 
effects. As shown in Ref. [ITJ, model-dependent acceptance effects can reduce the cross 
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section by as much as a factor of ~ 1/2 in the high bins. No uncertainties from the 
choice of scale or PDFs have been included either. It is reasonable to assume that these 
unaccounted for effects could help to ameliorate some of the tension between the measured 
and predicted value of {da / dMu) / cr . Nevertheless, the shape of differential cross section 
constrains the spin-2 model's parameter space. The question becomes, how large of an 
asymmetry can generated once this constraint is taken into account. 

V. CONCLUSIONS 

If it persists, the anomalously large top quark forward-backward asymmetry observed 
at the Tevatron is an indication of physics beyond the Standard Model. Since its ap- 
pearance, many models involving new scalar, as well as vector particles around the weak 
scale have been proposed to address the anomaly. However, it has proven to be hard 
to raise the theoretical prediction to the central value of the CDF measurement in the 
high mass bin. We have shown that there is parameter space in this model, consistent 
with various experimental constraints that could accommodate the CDF measurement of 
Ap^{Mtt > 450GeV) of 47.5%. The peculiar derivative coupling of a spin-2 particle to 
fermions naturally leads to strong sensitivity of the asymmetry to the ti invariant mass. 
As a result, the picture of the top asymmetry increasing with energy observed by CDF 
naturally fits in this framework. If the observed A'^p^ holds, it would be interesting to 
study the experimental bounds as well as the phenomenology of this model in more detail. 
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Appendix A: Form Factors 



The combinations of coefficients that appear in Eqs. (13), (14), (15) and (16) are 



^i — \9ut\ "T \ytu\ ~ \9ut\ " \9tu\ 
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C2 

C4 
C5 



ut\ ^ \ytu\ ^ \yut\ ^ \ytu\ ■> 

/j|2 I L|2\/| /j|2 I ^|2 
5'tju I + 1 9uu I ) \\9tt\ + I fi'tt I 

^««r+kir)^«(^a*+c.c., 



R |2 







iftir 






9ut\^ + \9tf) 9uu {9uT + {\9ut\^ + kf.r) ^i (^«)* + c.c, 
^^tf + kt^r) ^™ {9ttT + (kutf + kt^r) 9uu {9uT + c.c. 



The form factors in Eqs. (13), (14), (15) and (16) are: 



Fi(x, 


y) 


F2(x, 


y) 


Fz{x, 


y) 


Fi{x, 


y) 



F^ix, y) 



Ffs{x, y) 



Fj{x, 


y) 


Fs{x, 


y) 


Fg{x, 


y) 


Fw{x, 


y) 


Fuix, 


y) 


Fuix, 


y) 


Fisix, 


y) 


Fi4(x, 


y) 



{2y-2- x) (2y^ + Axy{A + x) - y'^{Q + 5s) - a;(8 + 8a; + x^)) , 

-y [{2y^ + Axy{A + x) - y^{Q + 5x) - x{8 + 8x + x^)) 

+ {2y~2-x) [y^ - 2xy + x{2 + x))] , 

y^ [y^ - 2xy + x{2 + x)] , 

18 (32 + 64x + A2x^ + lOx^ + x^) - 36y (64 + 82a; + 28^^ + 3x^) 

+ y'^ (2712 + 1764a; + 235a;2) - 2^^(492 + 109a;) + 73|/^ 

-2y [-3 (192 + 276x + 104^^ + 7a;^) + 2y (672 + 480x + 43x2) 

-1/2(720 + 109x)+44y3] , 

3 [3x^(2 + xf - 2xy (48 + 52x + llx^) + y'^ (256 + 332a; + 65x2) 

-4y^(72 + 19x) + 30y^] , 

-8y [-3x2(2 + x) + 2x2/(12 + 7x) - y^{2A + 19x) + 8y^] , 

162/2(x - ?/)2, 

6(1 + x)2(4 + x) - y (48 + 78x + 23x2) + y^{37 + 28x) - lly^ 

3x2(2 + x) + y (24 + 3x - 8x2) ^ y2^j^ - 8) - 2y^, 

Ay [x2 + y{l - 2x) + y^] , 

1 + lOx + 42x2 + 64x^ + 32x^ - 2y (3 + 28x + 82x2 + 64x^) 

+ 22/2 (9 + 68x + 96x2) _ ^^s^g _^ ^^x) + 32y^, 

2y (1 + 2x + 2x2 - 42/(1 + x) + 2y^) , 

(1 + 2x - 2y) (1 + 8x + 8x2 _ ^y^^^ ^ 4^^^ ^ g^2^ ^ 
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Fi5(a;, y) = 6(1 + xf (4 + 17a; + Ax^) - 2% (5 + 22x + 22^^ + 5x^) 

+ y'^ (253 + 596a; + 216a;^) - 2|/^(109 + 84a;) + A8y\ 
Fie{x, y) =y[3{8 + 33a; + 32a;2 + 8x^) - 2y (40 + 95a; + 36x^) + ^^(94 + 72a;) - 24y^] 
F,rix,y)=Ay\l + 2x-2y), 

Fis{x, y)=y [24 + 78a; + 51^^ + 4x^ - Ay (l8 + 28a; + 3a;^) + y^(61 + 12a;) - Ay^] , 
Fig{x, y) = 3x^ (6 + 11a; + Ax^) - Axy (l2 + 36a; + 17a;2) + y^ (24 + 191a; + 132a;2) 

-Ay^{20 + 27x) + 32y*, 
F2o{x, y) = 4y{x - y) (a;(3 + Ax) - y{l + 8a;) + Ay^) . 
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